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Abstract. Frequency-entangled photons can be readily produced using parametric 
^\ , down-conversion. We have recently shown how such entanglement could be 

manipulated and measured using electro-optic phase modulators and narrow-band 
frequency filters, thereby leading to two-photon interference patterns in the frequency 
domain. Here we introduce new theoretical and experimental developments showing 
f»^ ■ that this method is potentially a competitive platform for the realization of quantum 

communication protocols in standard telecommunication fibres. We derive a simple 

theoretical expression for the coincidence probabilities and use it to optimize a 

Bell inequality. Furthermore, we establish an equivalence between the entangled- 

k^( , photon scheme and a classical interference scheme. Our measurements of two-photon 

5_j ' interference in the frequency domain yield raw visibilities in excess of 99%. We use 



a 



our high quality setup to experimentally validate the theoretical predictions, and in 
particular we report a violation of the CH74 inequality by more than 18 standard 
deviations. 



PACS numbers: 42.65.Lm, 03.65.Ud, 03.67.Bg 
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1. Introduction 

Precision manipulation of entangled photons is highly desirable, both from the 
fundamental point of view of studying the ultimate limits of optics and from the point 
of view of applications such as quantum communication. Indeed, since Ekert's seminal 
work [1] , entangled photons have appeared to be a promising way to distribute quantum 
information. Using entangled photons could potentially allow the realization of key 
distribution protocols over distances greater than a few hundred kilometres [21 El 0] and 
security certification without a priori trust in the devices employed [5]. 

Most practical quantum key distribution methods based on entangled photons 
use time-bin [BJ [7] or polarization [SI [9] encoding. These have also been among 
the preferred methods for investigating the fundamental issue of quantum nonlocality 
[TUl [TTl [T2] . Manipulating entangled photons directly in the frequency domain is a 
relatively unexplored area. Previous work in this direction includes Hong-Ou-Mandel 
dip experiments p31 El [15], creation of entanglement in multiple degrees of freedom 
including frequency [161 E] and conversion from polarization to frequency entanglement 

In [19], we introduced the notion of frequency-bin entanglement that allows a simple 
description of experiments that manipulate entanglement in the frequency domain. 
We have shown how, using conventional methods of production (parametric down- 
conversion) and detection (avalanche photodiodes (APDs)), frequency-bin entangled 
photons at telecommunication wavelengths (about 1550 nm) could be manipulated in 
optical fibres using standard telecommunication components such as fibre Bragg gratings 
and electro-optic phase modulators (EOPMs) driven by radio-frequency (RF) signals. 

In this work, we improve on the work reported in [19] . We develop the theory behind 
this experiment as well as report experimental improvements. The results reported in 
[19] and in the present paper build on earlier experimental investigations on quantum 
communication using attenuated coherent states and frequency encoding [201 [21], which 
in particular aimed at applications in quantum key distribution. Complementary 
theoretical studies of the manipulation of photons in the frequency domain using EOPMs 
can be found in [22l [23] . 

The first part of the paper deals with the theoretical aspects of our experiment 
involving frequency-bin entangled photons and EOPMs. After describing the basic 
scheme, we show that it is possible to derive a considerably simpler expression for 
the joint probabilities than was reported in [19]. This shows that, whereas in most 
experiments involving entangled photons the interference pattern is given by sine 
functions, here it is given by Bessel functions. We then use this expression to derive the 
optimal RF settings for violation of a Bell-type inequality, the CH74 inequality. Finally, 
we adapt the correspondence between "prepare-and-measure" and entanglement-based 
schemes (often exploited in quantum key distribution) to our setup, showing that 
identical RF settings will give rise to identical interference patterns both in two-photon 
and in single-photon experiments. This identity is very useful experimentally as it makes 
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it possible to test the quality of the RF setup using a broadband white light source. A 
detailed comparison of the different experiments is reported in section [31 

At the experimental level, we have improved the experiment reported in [19] in 
a number of ways. We have developed an RF architecture based on off-the-shelf 
components that provides highly stable control of both the amplitude and phase 
of the RF signals used to drive the EOPMs. In addition, we have improved the 
stability of the pump laser wavelength and optimized the conversion efficiency of 
the periodically poled lithium niobate (PPLN) crystal. We now also use low-noise 
superconducting single-photon detectors (SSPDs). Altogether, these improvements 
allow us to report two-photon interference with raw and net visibilities of (99.17±0.11)% 
and (99.76 ± 0.11)%, respectively. These are comparable to the best results reported 
for two-photon interference at telecommunication wavelengths; see, e.g., [24]. By using 
the optimal settings mentioned above, we also report violation of the CH74 inequality 
by more than 18 standard deviations. Both these figures are considerably better than 
those reported in [19] . 

These results lay the groundwork for future experiments. In particular, they show 
that manipulating frequency-bin photon entanglement with EOPMs is a promising 
platform for the realization of quantum communication protocols at telecommunication 
wavelengths. 

2. Theoretical analysis 

2.1. Manipulating and measuring frequency- entangled photons using electro-optic phase 
modulators and narrow-band frequency filters 

The experimental scheme we study in this paper is depicted in figure [T^b), see section 

3. It consists of a continuous narrow-band laser at frequency 2uq pumping a parametric 
down-converter. After removal of the pump beam, the signal and idler photons are 
separated. They pass through EOPMs and through narrow-band frequency filters, and 
are finally detected. We now give a detailed theoretical description of this experiment. 

The parametric down-converter produces photon pairs in a frequency-entangled 
state of the form 

r+oo 

|tf>=/ duf{u)\u + u)\uj -uj), (1) 



where f(co) = r(ui)e 1 ^ U1 ^ is a complex function of uj that characterizes the bandwidth 
of the signal and idler photons. We have neglected in this expression the linewidth of 
the pump laser -- whose effect is to make the pump (angular) frequency 2o; slightly 
uncertain. 

If one measures the frequencies of Alice and Bob's photon state, one finds perfect 
correlations: if Alice obtains uq + u, Bob obtains uq — u. In practice the frequency can 
only be measured with precision Qp given by the width of the frequency filters used. 
This leads to the notion of frequency bin: all photons whose frequencies are contained 
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in an interval [co F — y,wp + Qp] are grouped into a single frequency bin centred on 
frequency wp. 

When a single photon of frequency u passes through an EOPM , driven at the RF 
f2 RF (with Q RF > Q F ) with adjustable amplitude c and phase 7, it undergoes the unitary 
transformation 

|w) H> U(c, 7) |w) = J2 Un ^ 7) l w + ^Rf) , (2) 

where 

f/n(c,7) = ^(c)e in ^/ 2 ) (3) 

and J n is the nth-order Bessel function of the first kind. 

Since the action of EOPMs on this state can only change the frequencies by integer 
multiples of fl RF , it is convenient to rewrite the state as 

/+n RF /2 
du ^ f( n ^RF + w) \u + ntt RF + to) \u - ntt RF - to) . (4) 

The motivation for re-expressing ([1]) in this form is that the EOPMs will cause 
contributions from different values of the index n to interfere, while contributions from 
different values of the offset parameter u will add probabilistically. Indeed, with a 
sufficiently precise measurement of the frequencies of the photons exiting the EOPMs, 
we could determine a specific value for uj and, in retrospect, claim the initial entangled 
state was 

Y$>) = y]f(nfl-RF + u>)\wo + nQ RF + u)\uj - nQ RF - ui) . (5) 

nGZ 

We can further simplify notation by noting that the actual value of u in (jHJ) is of 
no importance. We therefore drop the parameter u from (jSJ) and adopt the discretized 
version 

i^hE^i-™)' ( 6 ) 



iiRP I iZ 



RF 



where \n) denotes a photon with a frequency u + n£l RF + u for some u G [ — I s -, 4 
and we denote f n = r n ^ n = f(nQ RF + u>), Vn. 

We will also make the hypothesis that /„ varies slowly with n, which is justified 
if Q RF is very small compared to the frequency range over which / varies. In our 
experiments the bandwidth of the photon-pair source (the scale over which / changes) 
is approximately 5 THz, while f2 RF = 25 GHz. This allows us to identify f n & f n+p for 
small values of p, say — 5 < p < +5 (see next paragraph). 

In the experiment schematized in figure [T^b), each photon is separately modulated 
with respective parameters a, a and b,(3. According to (j2j), the state (J6]) is transformed 
to 

\*)^U(a,a)®U(b,P)\V)= J2 fnc d (a,a;b,(3)\n)\-n + d), (7) 

n,d£Z 
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where c d (a,a;b, (3) = ^ z U p (a,a)U d - p (b, (3), and where we use the assumption that 
f n ± p « f n . This is reasonable since \U p (c, 7)! decreases rapidly with p for accessible 
RF amplitudes: the values of p for which \U P \ is large are limited to approximately 

pe [-5, +5]. 

The joint probability of Alice detecting a photon in the frequency bin n on which 
frequency filter Fa is aligned and Bob detecting a photon in bin — n + d on which filter 
Fb is aligned is given by 

P d (a,a;b,l3;n) = \(n\(-n + d\q)\ 2 = \f n \ 2 \c d (a,a;b,(3)\ 2 . (8) 

At this stage, we note that the series giving the coincidence probability can be 
summed by using the Graf addition formula [25] (which we rederive in Appendix A in 
our adopted notation). We have 

J2 U P (a, a)U d . p (b, (3) = U d (C, T) , (9) 



h 2 J- %/irncL _ /^l 1 / 2 Qnrl tan F - 

a cos a+b cos B ' 



where C = [a 2 + b 2 + 2a6cos(a - (3)} 1 / 2 and tanT = ggg^| , and we may always 
take C to be positive. In terms of these parameters, 

P d (a,a;b,(3;n) = \f n \ 2 J d (C) 2 . (10) 

Thus manipulating frequency-entangled photons with EOPMs gives rise to Bessel-type 
interference patterns, rather than the usual sine and cosine interference patterns in 
optics when only two modes are present. 

Note that ( ITOl) implies the normalization 

J2 Pd(a, a; 6, /?; n) = \f n \ 2 ]T J d (C) 2 = \f n \ 2 (11) 

d d 

required by conservation of probability. 

Note also that with modulation turned off the photons do not change frequency 

and we have P d= o(a = b = 0; n) = \f n \ 2 and Pd^o(a = b = 0; n) = 0, as expected. 

Equation ( fTUl) shows that in the experiment schematized in figure [H(b), the 

(2) 
coincidence rate N d for frequency bins n and — n + d will be given by 

N d '(a, a; b,/3; n) 

= J d {[a 2 + b 2 + 2abcos{a - (3)} 1 ' 2 ) 2 x N ( d %(a = b = 0;n) , (12) 

(2) 
where N d J (a = b = 0; n) is the coincidence rate for frequency bins n and —n when the 

modulation is off. 

2.2. Bell inequality optimisation 

We now show that the correlations ( TTOl allow the violation of a Clauser-Horne (CH74) 
inequality [26], specifically the violation of S < 2, where (see [T9] ) 

(n\ /o\ /o\ 

S = [N K d J (a , a ; b , (3 ; n) + A^J (a , a ; 61, 0i; n) + N^J (a 1 , «i; b , P ; n) 
- N^%(a u au b u /3 i; n)]/^2 (o = 6 = 0; n) . (13) 
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The summation based on the Graf addition formula allows a relatively 
straightforward determination of the optimal parameters for violating the CH74 
inequality. To this end, we first substitute ( TTOj) in order to rewrite the CH74 expression 

as 



S = Jo(Coo) 2 + Jo(Coi) 2 + Jo(Cio) 2 - Jo(Cn) 2 , (14) 



where 



dj = [a 2 + b) + 2a i b j cos(ai - /^)] 1/2 (15) 

and i,j e {0, 1}. 

The parameters Cy obey constraints imposed by the form of (TT5|) . To see this, we 
introduce the vectors 

a,i = (ai cos Oj, Oj sin Oj) , 

bj = — (bj cos (3j, bjsia/3j). (16) 

In terms of these, Cy = |a, — 6j|. We may therefore identify the Cy with the lengths of 
the sides of a quadrilateral defined by the vertex vectors cto, &o 5 a i an d bi. This implies 
that each of the four CV,- is bounded by the sum of the other three. For example, 

Cn = |oi - 6i| 

= |oi - 6 + &o - «o + «o - bi| 

< |oi - 6o| + l^o - a o| + l a o - bi\ 

= C W + C 00 + C 01 . (17) 

In this way, we reduce the eight-parameter optimization of (1131) to a four-parameter 
optimization with constraints. There are two possibilities: either the optimum will lie 
within the parameter domain, or it will lie along one of the boundaries. We quickly 
rule out the former possibility: in this case, local extrema are found for parameters 
Cij which are local extrema of J 2 . Since at the second greatest extremum we have 
Jo(x) 2 ~ 0.162 (for x ~ 3.832), no combination of four positive extrema of J 2 satisfying 
the strict inequalities will lead to a violation of the CH74 expression. 

A Bell inequality violating optimum, if one exists, must therefore lie along one of 



the boundaries. In Appendix B we show that the global optimum of (IT4"|) lies along the 
constraint C 00 = C m = C 10 = Cn/3, systematically ruling out any other possibility. 
The optimum of Coo ~ 0.550, for which S ~ 2.389, corresponds to the RF parameters 

(ao,a ) = (0.275,0) =(&o,A>), 

(oi, ai) = (0.825, + tt) = (&i, A) . (18) 

Reaching the optimal value requires the use of variable - - but small - - modulation 
amplitudes and precise phase adjustment. 

These optimal parameters should be contrasted with the analysis and experiment 
reported in [19] where the CH74 inequality was violated in a configuration where the 
amplitudes a = Oi = b = b\ were all equal, and only the phases a , «i, (3 , fix varied. In 
this configuration, the largest CH74 violation tends to the above, but requires arbitrarily 
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large RF amplitudes. But when the amplitudes are large, the two-photon interference 
pattern is much more sensitive to small errors in the RF amplitudes and phases. The 
above — optimal — value of S is attained for rather small values of the RF amplitudes, 
which makes the experiment much more robust. 

For comparison, we note that the maximum value for the expression ( fT31 attainable 
by quantum theory for systems of dimension 2 is 2.414, which is quite close to 
the maximum value of 2.389 attainable using EOPMs on frequency-bin entangled 
photons. However, as the frequency-bin entangled photons belong to a Hilbert space 
of dimension greater than 2, it may be that the maximum value attainable by some 
local measurements on the state exceeds this value. We do not know whether this is 
the case. The algebraic maximum for this expression is 3, which cannot be exceeded by 
any measurement. 

2.3. Equivalence of two-photon and single-photon interference schemes 

There is a mathematical correspondence between correlation experiments on maximally 
entangled states and prepare-and-measure schemes, based on the identity (i\(j\Ua <8> 
f/ B |$+) = (i\U A U£ \j)/y/d where |$+) = Ya=i \i)\i)/Vd- Indeed, the first term in the 
equality can be interpreted as a measurement on the entangled state |<3? + ) in which 
Alice projects onto (i\U& and Bob projects onto (J\Ub, whereas the second term can be 
interpreted as the preparation by Bob of the state U^\j) and the subsequent projection 
by Alice onto the state (i\Ua_- 

This theoretical correspondence is well established in the context of quantum key 
distribution, where it is used to demonstrate the equivalence between prepare-and- 
measure schemes and entanglement-based schemes. However, in general, these two 
schemes correspond to different experiments because implementing U£ is physically 
different from implementing U&. This is the case, for example, in experiments involving 
time bins. In our case however, where the transformations £/a,b are realized by EOPMs, 
Ub = Ub, in the sense that we have the identity (—n\U(a,a)\p) = (—p\U(a,a)\n). 
The above mathematical identity thus translates to a physical correspondence between 
transition amplitudes in two-photon and single-photon experiments, with all the RF 
parameters (amplitudes and phases) kept unchanged. 

In practice, this corresponds to the equivalence between the scheme depicted 
in figure [H(b) in which entangled photons are manipulated by EOPMs and that of 
figure QJa) in which photons belonging to a particular frequency bin are selected by 
filter Fa and subsequently modulated with parameters (a, a) and (b,(3). Specifically, 
the amplitude of detecting photons in frequency bins — n and n + d in experiment QJb) 
is proportional to the amplitude of detecting a photon in frequency bin n + d given that 
it was prepared in bin n in experiment Ufa): 

(-n\(n + d\U(a, a) ® U(b, /3)|tf) 
= Y,Mn + d\U(b,P)\-p)(-n\U{a,a)\p) 
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^ f n J2( n + d P(b^)\-p)(-p\U(a,a)\n) 

p 
= f n (n + d\U(b,P)U(a,a)\n), (19) 

where in line 3 we have invoked the same assumption used to derive dTJ) , namely that /„ 
is approximately constant over the experimentally accessible range of frequencies which 
interfere, and we used the completeness relation ^2 \p)(p\ = 1. 

Thus, in the experiment of figure [T^a), if a single photon is initially in bin n, the 
probability that it is detected in bin n + d has the same functional dependence as 
( TTOl) in the two-photon case. Experimentally, the quantity measured is the photon rate 
AQ , where d denotes the separation between the initial bin n and the final bin n + d. 
Therefore, we have 

N d (a, a; 6, /3; n) 

= J d ([a 2 + b 2 + 2abcos(a - /3)} 1/2 ) 2 x N^ (a = b = 0;n) , (20) 

where N d J (a = b = 0; n) is the photon rate for frequency bin n when the modulation 
is off. 

Moreover, if bin n initially contains a large number of photons, then in the 
experiment of figure [T]^a) the optical power N d measured in bin n + d, if the initial 
light beam is prepared in bin n, is given by 

N d c (a,a; b, f3;n) 

= J d ([a 2 + b 2 + 2abcos(a - (5)} l/2 ) 2 x N^ s \a = b = 0; n) , (21) 

where AQ =0 (a = b = 0; n.) is the optical power for frequency bin n when the modulation 
is off. 

In addition to the correspondence to the entanglement-based scheme, experiments 
based on figure HJa) are of interest in themselves: when used with a single-photon 
source, this scheme allows the realization of quantum communication protocols such as 
quantum key distribution, see [2T] . 



2.4- A note on phase accumulation during propagation 

If a photon of frequency uj propagates a distance L, its state is transformed according 

to 

where we have developed the wave number k in series in u — ojq. 

When considering frequency coding of information, the zeroth-order term e 1/3 ° L is 
an overall phase, and has no physical influence. The first-order term q 1 Pi( u - u >o)l can 
be absorbed into the phase 7 of the RF field applied to the EOPM, see ([3]). In all 
calculations as well as the analysis of the experiments, we absorb the phases e 1 ^ u '~ Wo ' L 
into the phases of the RF fields. 

The fact that phases accumulated during propagation can be absorbed into the 
phases of the RF fields underlies the inherent high stability of experiments using 
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frequency bins. Indeed, stability of our experimental setup requires that /^if^RpL <C 1 
(since the only frequencies that interfere are those separated by small multiples of 
Hrf). This should be contrasted with interference experiments in the spatial domain, 
where one is sensitive to the phase (3 L. Approximate equality of the phase and group 
velocities implies that /3 ~ PiuJq. Since Q^p/uo — 10~ 4 , our experiments are less 
sensitive to changes in fibre lengths by a factor of roughly 10 -4 . This implies that 
in laboratory experiments, no stabilization is required. In field experiments, however, 
where propagation distances are tens of kilometres, the local RF oscillators must be 
synchronized, as was done for instance in [27]. 

Note that the higher-order terms in (1221) due to frequency dispersion cannot be 
absorbed in the phase 7 of the RF field. In this work, we neglect dispersion effects, but 
it may degrade the quality of interference in certain cases — i.e. dispersion management 
may be needed for long-distance and broad-spectrum applications. 

3. Experimental setups 

In this section we describe the experimental implementations of the setups already 
introduced and illustrated in figures UJ^a) and (b), and point out the critical requirements 
for high-visibility experiments. 

3.1. Classical optics and one-photon experiments 

We present here in total three versions of the experiment schematized in figure H](a), 
differing mainly in the optical source S employed. In the first version, S is a classical 
source consisting of a coherent polarized narrow-band laser. In this case neither filtration 
nor polarization management is required, since all the photons emitted already belong 
to a given frequency bin, with their polarization aligned with the active axes of the 
modulators. In the second version, the source, also classical, consists of a non-coherent 
non-polarized broadband light source. In this case, the light beam must first pass 
through the filter Fa and a polarizer before being sent through the setup. This 
corresponds more strictly to the entangled-photon case, where signal and idler photons 
have a broad spectrum. In the third version (the one-photon experiment), an ideal 
broadband single-photon source is approximated by attenuating the broadband light 
source until it contains, on average, much fewer than one photon in each frequency bin 
within each detection time window. We report the results for all three sources. 

In the experiment, filter Fa is a home-made fibre Bragg grating used in reflection 
with a circulator. It selects photons belonging to the frequency bin [cjf — y,wp + -if\ 
centred on u-p = coq and restricted to a 3 dB width -^ w 3 GHz. Photons are 
consecutively guided through EOPM A and EOPM B (EOSPACE) driven by an RF 
signal at ^^- = 25 GHz with adjustable amplitude and phase a, a and 6, /3, respectively 
(controlled by the setup depicted in figure |2]). At 12.5 GHz, the isolation of Fa is better 
than 30 dB, ensuring that frequency bins are well isolated from each other. Photons 
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Figure 1. Schematic representation of the experimental setups, (a) Classical optics 
and single-photon experiments. The source S can be a narrow-band laser, a broadband 
incoherent source or a broadband incoherent source attenuated to reach the single- 
photon regime. The detector D is either a classical photodiode or, in the single-photon 
regime, an APD operated in the Geiger mode, (b) Two-photon experiment. The source 
consists of a PUMP laser (whose wavelength Ap is stabilized) that enters a PPLN 
waveguide (whose temperature Tp is stabilized) and is then removed by a drop filter 
F. Detectors Da and Db are SSPDs. In both panels, band-pass filters Fa and Fp select 
photons belonging to a specific frequency bin. The phase of the light is modulated by 
EOPMa and EOPMb , driven by a signal at frequency f^ap with respective parameters 
a, a and b, (3, determined by the RF system depicted in figure [2] 



selected by filter Fb (which has the same characteristics as Fa) are detected. Where 
classical sources are used, the optical power j\K class ) is measured by a photodetector D. 
In the single-photon regime, the single-photon rate N^ is measured with an APD (id 
Quantique) operated in gated mode (repetition rate = 100 kHz, detection window size 
= 2.5 ns, efficiency ~ 10% and dark count rate fa 3 x 10~ 6 ns -1 ). 



3.2. Two-photon experiment 

In the experiment schematized in figure QJb), a narrow-band (2 MHz bandwidth) 
continuous pump laser (Sacher Lasertechnik) stabilized at a wavelength Ap emits 3 mW 
of power into a 3 cm long PPLN waveguide (HC Photonics). Phase matching is achieved 
by controlling the waveguide's temperature T P . It is then possible to efficiently generate 
frequency-entangled photon pairs centred on a specific frequency, corresponding in our 
case to the wavelength 2Ap = 2tic/uj = 1547.743 nm. 

At the output of the PPLN, a drop filter F rejects the pump with more than 60 dB 
isolation. The pairs are separated (with 50% probability) with a 3 dB coupler, such that 
each photon in every entangled pair is sent through an independent EOPM, EOPMa 
or EOPMb- Polarization- maintaining fibre components ensure that the photons' 
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polarizations are aligned with the active axes of the modulators. Filters Fa and Fb 
select photons belonging to the frequency bin centred on uj . Photons are finally detected 
by SSPDs Da and Db (Scontel) cooled to 1.7 K and operated in the continuous mode 
(efficiency ps 5%; dark count rate ps 30 Hz). 

A time-to-digital converter performs a coincidence measurement. Specifically, it 
registers the arrival times £a and £b of photons A and B and records the number 
of coincident detection events as a function of £b ~ ^a- The coincidence rate N^ is 
extracted from the histogram by summing all contributions for which the difference in 
arrival times £b — ^a is contained in a given time interval of width 0.6 ns. 

3.3. Requirements for high-visibility experiments 

In order to produce high-visibility Bessel interference patterns, our architecture must 
fulfil some critical requirements, particularly in the two-photon case. 

First, in order to obtain precise reproducible interference in the frequency domain, 
high-resolution phase and amplitude control of the RF signals must be achieved. We use 
the RF architecture presented in figure El This architecture is based on RF translation 
and uses cheap off-the-shelf RF components. It could therefore be adapted for field 
experiments in which Alice and Bob are separated by a large distance. 



ALICE 




Figure 2. The RF system used to drive EOPMa and EOPMb in all experiments. A 
2.5 GHz signal of an RF generator G is separated by power splitters PSq, PSa and PSb- 
The phase of both Alice's and Bob's signals is controlled by I&Q modulation. Using 
frequency multipliers x4 and x2, signal multipliers Ma and Mb and amplifiers Aa 
and Ab, strong 25 GHz signals are obtained. Their amplitude is controlled by variable 
attenuators. Accurate control of the parameters a, a and 6, /3 is thus obtained. 
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In this RF architecture, an RF generator outputs a 2.5 GHz, 10 dBm signal. Its 
power is equally shared by splitter PSg between Alice and Bob. The two parts of the 
setup are equivalent and made independent with isolators. Alice's and Bob's signals 
are split by PSa and PSb, respectively. A first part of the signal is subjected to I&Q 
modulation: a given combination of attenuation allows the precise imposition of an RF 
phase. A digital-to-analogue converter controls with high precision the I&Q voltages 
Ia-iQa and Ib,Qb, allowing fine selection of (I,Q) pairs. In this way, phases a and 
ft can be shifted from to In while the amplitude remains constant. Since the phase 
shift introduced at 2.5 GHz is altered by frequency translation (see below), a mapping 
is performed at 25 GHz in order to identify the correct (/, Q) couples. This permits 
distortion correction, resulting in a highly accurate control of the phase of the RF signal 
with a precision better than 10™ 2 rad. The second part of the signal is sent through a 
multiplier and the 10 GHz signal obtained is multiplied, with fixed phase, by the phase- 
shifted signal to obtain a 12.5 GHz signal with adjustable phase. Frequency filters 
control the purity of the signal. The signal is sent through an amplifier and a frequency 
doubler, which amplifies the signal and converts it into 25 GHz. Finally, the signal 
is applied to EOPMa or EOPMb with an adjustable amplitude a = nV^ 1 R 1 ^ 2 P A or 
b = nV~ 1 R 1 / 2 P B , where the half-wave voltage V^ ~ 2.9 V, the internal resistance 
R = 50 Q and power Pa or Pb is mechanically adjusted with a variable attenuator. Ten 
per cent of the signal is sent to an RF power meter allowing an adjustment of Pa and 
Pb with 10~ 2 dB precision. 

Second, in order to achieve a high and stable SNR in the two-photon case, fine 
control of the pump laser wavelength is compulsory. Each down-converted photon 
must belong to a frequency bin whose width is set by the 3 GHz band-pass filters 
Fa and Fb- The pump wavelength must thus be accurately set at the centre of a 
frequency bin. To this end, 10% of the power of the pump laser is sent to a wavelength 
meter (Exfo) that generates an electric signal proportional to the difference between the 
measured wavelength and a reference. A proportional-integral-derivative loop generates 
an error signal which feeds the piezoelectric transducer of the external cavity of the laser, 
stabilizing the wavelength at Ap = (773.8715 ± 0.0002) nm. The degeneracy frequency 
ouo is thus controlled with a precision of around ±0.04 GHz. 

Third, in order to have a high and stable pair production rate, one must also 
optimize the parametric down-conversion of the PPLN source. To achieve a fine 
conversion wavelength control and good efficiency, the PPLN crystal is seated on a 
Peltier cell for accurate temperature stability. Precise measurements with a tunable 
C-band laser in the second-harmonic regime show that after several hours, a stable 
optimized conversion wavelength is reached with about 1 pm precision. 

Achievement of the latter two conditions allows us to detect (using low- noise SSPDs) 
coincidences at a rate of approximately 20 Hz and with a SNR of approximately 2 x 10 3 
(without RF signal applied to the EOPMs). 

Together with the accurate control of the RF parameters, these improvements lead 
to raw and net visibilities of (99.17 ± 0.11)% and (99.76 ± 0.11)%, respectively. This 
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should be compared to the 98% net visibilities reported in 

Note that compared to our previous work [19] . control of the RF phases is now 
automated and the accuracy much improved. The detectors in all experiments can be 
interfaced to a computer, enabling both selection of RF phases a and (3 and measurement 
data acquisition. Rates N^ 2 \ N^ and ]\[( cla - ss ) are thus measured during automatically 
adjustable times in automatically adjustable RF configurations. Note that replacing 
mechanically variable attenuators by remotely controlled components would allow full 
automatization of the system. Note also that for long-distance experiments it is possible 
to use two synchronized RF generators instead of one, as shown in [27] . 

4. Experimental results 

4-1. Equivalence of interference schemes 

We first evaluate the rates N^ 2 \ N^ and jv( class ). As discussed in section 12.3] 
these should have identical dependence on the RF parameters, see (TT2l) . ( |20l and 
( 12T1) . To experimentally test this equivalence, we chose d = 0, and the parameter 
C — [a 2 + b 2 + 2abcos(a — /3)] 1//2 is varied by scanning one of the phases a or (3 with 
a = b fixed. This procedure permits easy evaluation of the interference visibility, the 
value of which is critical for the performance of the system. 
We define the raw and net visibilities as follows: 

t r _ -* 'max J- >min /OQ"! 

Kraw ~ n +N ■ ' l ' 

1 ' max T^ J ' mm 
-,7- _ I,-' ''max -'''noise J l^'min -'''noise/ ( r )A\ 

= JN - N ■ ) + (N ■ - N ■ ) ' [ ' 

V iv max iv noisej T^ V/'min iv noisej 

where N denotes either N^, N^ or J\[( class )^ and A^ no ; se in (124]) represents the noise due 
to detector imperfections (e.g. dark counts) which can be measured independently. 

The maximal rate iV max is obtained in principle for C* = 0, such that Jo(C*) 2 = 1. 
For a = b, this is achieved with a phase difference a — (3 = n. The minimal rate A^ min 
is obtained for any of the positive roots {C*} of J , for which J (C*) 2 = 0, Vi. The 
visibility therefore attains a theoretical maximum value of 1. If a = b, the first root 
C* is attainable at sufficiently high RF powers (a = b > 1.2) with the phase difference 
a -(3 = arccos(C 1 * 2 /2a 2 - 1). 

Our results are summarized in figure [3] and in table [U Figure [3] is a plot of 
Jd{[a 2 + b 2 + 2abcos(a — /3)] 1 / 2 ) 2 as a function of a — (3 for d = and a = b w 2.25. 
The experimentally derived rates N^, N^ and J\[( class ) - - normalized by iV- 2 - 1 , iV^ 1 ) 
and j /\T ( ' class * ) with modulation turned off — are superposed to the theoretical curve. The 
close agreement of the experimental data taken with the different experimental schemes 
validates the predictions of ( IT2l) . ( 1201) and (12TT) and thereby demonstrates the equivalence 
of the two-photon, one-photon and classical optics experiments of figure [TJ 

Note that, as mentioned in section [31 an uncontrolled phase shift arises due 
to propagation times in the optical fibres. The experimental results are therefore 
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Figure 3. Bessel interference pattern in the frequency domain. The interference 
described by (|T2j) . (|20|) and (|2Tj) is plotted for d = and a = 6 f=a 2.25 as a function 
of the phase difference a — j3. The raw (with detector noise) and net (detector noise 
subtracted) data are shown, respectively, on the top and bottom panels. On both plots, 
the theoretical curve is included in black, the green data points are the normalized light 
intensity JV( class ) (obtained using the tunable laser source and the setup of figure [TJa)), 
the blue data points correspond to the normalized photon rate N*- 1 ' , i.e. the single- 
photon interference pattern (obtained using the broadband source, attenuated to the 
single-photon regime, and the setup of figure [TJ a)), and the red data points correspond 
to the normalized coincidence rate N^ 2 \ i.e. the two-photon interference pattern 
(obtained using the PPLN source and the setup of figure QXb)). The error bars are 
statistical. Note that the classical measurements 7V( class ) are made with a relatively 
noisy photodiode, and the raw visibility for the green data points in the top panel 
is therefore less than that reported in table [TJ where a low-noise photodiode was 
used. Note also that when detector noise is subtracted (bottom panel) the different 
experimental data points superpose exactly, demonstrating the equivalence of the 
interference schemes, but deviate slightly from the theoretical curve. We attribute 
this to small errors in the calibration of I&Q parameters, so that the actual phases a 
and /3 deviate slightly from their theoretical value. 
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horizontally shifted in order to obtain the best possible agreement with the theoretical 
curve. This is the only parameter that is adjusted to fit the data. 

Table [D gives the values of the visibilities V Taw and V net extracted from curves such 
as those reported in figure |3j In the classical case a very low noise detector was used to 
measure maximal and minimal rates, which is why V raw ~ V^t- Both the laser source 
and the broadband source were used. In the single-photon case an APD with a relatively 
high dark count rate was used, which explains the low value of the raw visibility. In the 
entangled-photon case, the high value of V raw is due to the quality of the SSPDs, which 
are much less noisy than APDs. Maximal and minimal rates were measured for several 
minutes in order to obtain good statistical precision on the visibilities. 

Table 1. Experimentally measured visibilities. 

Classical expt Classical expt One-photon expt Two-photon expt 

(laser source) (broadband source) 

Kaw (99.79 ± 0.01)% (99.41 ±0.12)% (87.25 ± 0.38)% (99.17 ± 0.11)% 

Ket (99.79 ±0.01)% (99.41 ±0.12)% (99.27 ± 0.43)% (99.76 ± 0.11)% 

The values of V Taw depend strongly on the noise inherent in the detectors and thus 
vary greatly across the different cases. By contrast, the values of V ne t are all almost 
equal and notably high. This agreement once again confirms the equivalence between 
the different experiments. It also allows us to separate the contributions to visibility 
degradation due to the experimental setup from those inherent to the detectors (detector 
noise and dark counts) and sources (in particular noise due to multiple photon events) 
used. The main contributions from the setup itself are imperfect frequency-bin isolation, 
imperfect polarization management and imperfect control of the RF parameters. 

4-2. Bell inequality violation 

As a second main result, we evaluate the experimental violation of the CH74 inequality 
introduced in section [221 Scans of both phases a and (3 at the amplitudes (ao,bo), 
(a , bi), (oi, b ), (a l5 6i) given in (fi~8l) enable precise selections of the phases a , (3 , a±, fi\ 
given in ( TT8l) for which the violation will be largest. Measurements are then realized for 
phases and amplitudes optimizing the violation. 

The experimental results, listed in table |2j agree with the theoretical predictions 
up to the statistical errors and imply that the CH74 inequality is violated by more than 
18 standard deviations. 

We should, however, mention that this result does not provide a decisive test of 
local causality as we have not closed either the detection or locality loopholes, and 
because the CH74 inequality, as we have applied it, requires additional assumptions (see 
the discussion in JT9J [26] ) . Nevertheless, these results show that our present approach 
allows the study of quantum correlations of frequency-entangled photons, and could in 
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Table 2. Bell inequality violation results. AiBj, i,j £ {0,1} is the notational 
shorthand for (eij, &i;bj, /3j). The column "Theory" gives the values that should be 
obtained at the point ([T8j) for which the CH74 violation is maximal. Note that the 
experimental values agree with theoretical predictions up to statistical errors. 





Theory 


Expt Nil 


Expt N™ 






(with noise) 


(noise subtracted) 


A B 


0.857 


0.862 ±0.006 


0.861 ±0.006 


A B 1 


0.857 


0.863 ±0.006 


0.862 ±0.006 


A 1 B Q 


0.857 


0.854 ±0.006 


0.853 ±0.006 


A 1 B 1 


0.182 


0.190 ±0.003 


0.186 ±0.003 


S 


2.389 


2.389 ±0.021 


2.391 ±0.021 



principle, i.e. if experimental imperfections (mainly losses and detector inefficiencies) 
were small, be adapted to permit a decisive test of local causality. 

5. Conclusion 

In summary, building on our earlier work [19] . we have further developed the theory 
underlying the manipulation of frequency-bin entangled photons with EOPMs, and 
demonstrated that this could be reliably realized experimentally. At the theoretical 
level, our main results are an analytic simplification of the expression giving the 
coincidence probabilities, and the demonstration of the equivalence of prepare-and- 
measure schemes with two-photon schemes. In our experiment, frequency-entangled 
photons are produced by parametric down-conversion with a frequency- stabilized pump 
laser and a temperature-stabilized PPLN waveguide. They are made to interfere in the 
frequency domain through EOPMs whose driving signals are controlled by a dedicated 
RF architecture built entirely from off-the-shelf components. Interference patterns are 
detected with SSPDs via optical frequency filtering. The interference patterns are 
accurately controlled, exhibit high visibilities and allow the violation of a Bell inequality 
by 18 standard deviations. 

The strengths of our method are 

• the use of optical, electro-optic and RF components that are commercially available 
and allow easy interconnection and remote control; 

• the use of optical components that allow good polarization management, frequency- 
bin isolation and stability; 

• the use of an RF system that allows stability, independence, easy calibration and 
precise adjustment of parameters; 

• overall reproducibility and robustness allowing day-long experiments with no 
measurable drift or decrease in performance. 

Our system can be easily adapted to field experiments, as it is robust and allows full 
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automatization and long-distance synchronization. We are therefore confident that 
we will be able to extend our results in the near future to long-distance quantum 
communication experiments and to perform tasks such as long-distance Bell inequality 
violation or quantum key distribution. These future experimental works would be 
supported by the theoretical analysis introduced here. 
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Appendix A. The Graf addition formula 

The matrix elements U n (c, 7) = J n (c)e m ( 7_7r / 2 ) are the Fourier components of the 
periodic function e _lccos ^ _7 \ so we have 

e - icC os(0- 7 ) = J2 U n (c, 7 )e- in ^ (A.l) 

n 

(the Jacobi- Anger expansion). Since adding two sine waves with the same period yields 
another sine wave, we may write 

e -iacos(0-a) e -ifecos(0-/3) _ Q -iC cos{4>-T) /^ 2) 

We extract an expression for C and V by applying the identity cos(0 — 7) = cos cos 7 + 
sin sin 7 and comparing terms in sin (p and cos 0: 

a cos a + b cos (3 = C cos Y , (A. 3) 

a sin a + b sin (3 =Csinr, (A. 4) 

from which we obtain 

C 2 = a 2 + b 2 + 2ab cos(a - 0) (A.5) 

and 

asina + 6sin/3 

tan T = — . (A.6 

a cos a + b cos p 



Note that we may impose C > 0, which fully determines I\ Equations ( 1A.1I) and ( 1A.2I) 
and the convolution theorem then imply 

J2 U n (a, a)U d ^ n (b, (3) = U d (C, T) . (A.7) 
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Intuitively, ( 1A.2J) and ( 1A.7J) express that two phase modulators used in series, driven 
by sinusoidal RF signals of the same frequency, have the same action as a single phase 
modulator. Equation ( 1A.7J) is a slight generalisation of the Graf addition formula as 



given in [25], expressed in the notation we have adopted. 

Appendix B. Optimization of CH74 

We begin with the boundary Cn = C 00 + C i + C 10 . Considering the partial derivatives 
of the expression 

S = J (C 00 ) + Jo(Coi) + ^o(Cio) — </o(Coo + C*oi + C*io) > (B-l) 

we find that the optimum must satisfy 

( J 2 )'(C 00 ) = ( J 2 )'(C 01 ) = (J 2 )'(Cio) = (</ 2 )'(Ox, + On + C w ) , (B.2) 
where («/ 2 )' * s the derivative of J 2 . Due to the symmetry in ( IB. II) . we may impose 
do < Coi < Cio without loss of generality. Since J 2 < 1, a violation of S < 2 will 
clearly require Coo < £2/3 ~ 0.878, such that Jo(Coo) 2 > 2/3, and Coi < £1/2 ~ 1.126, 
where Jo(Coi) 2 > 1/2. Because there is no overlap in the images (Jq)'{[0,X2/3\) 
and (</ 2 )' 0^2/3, ^1/2]), ( IB. 21) further imposes Coi < £2/3- Finally, (IB. 21) together 
with injectivity of (Jq 2 )' on the domain [0,22/3] impose Coo = Coi, so we require 
Coo — Coi < £2/3- 

We first assume Coo = Coi = Cio < £2/3 and Cn = 3Coo, reducing the problem to 
a one-parameter optimization of 

S = 3J (Coo) 2 -Jo(3Coo) 2 . (B.3) 

Numerically, we find an optimal value of S'opthnized ~ 2.389 for C o ~ 0.550. 

We now show that this is the global optimum. The alternative would be to take 
some Cio > x i- in this case, it will not be possible to satisfy (IB. 21) unless Ci > 1.291, 
and we have no chance of surpassing the optimum given above unless Jo(Cio) 2 > 0.389, 
which translates to Cio ^ 1.293. Via ( IB. 21) . this in turn imposes Coo = Coi > 0.876, for 
which 2J (Coo) 2 < 1.634. However, J (Ci ) 2 < 0.390 for 1.291 < C w < 1.293, which 
guarantees we will not surpass the optimum of ( IB. 31) . 

Finally, we check that there is not a better optimum along one of the three other 
possible boundaries. Using the symmetry of (TT4l) . without loss of generality we consider 
the boundary Cio = Coo + Coi + Cn. In this case, and assuming that the function J 2 
is symmetric, the analogue of (IB.ll) is 

S = Jo(Coo) 2 + Jo(Coi) 2 + Jo(Cio) 2 - Jo(C 00 + On - Cio) 2 (B.4) 

and we find that an optimum along this boundary must satisfy 

(Jo 2 )' (C 0Q ) = ( J 2 )'(Coi) = (J 2 )'(-Cio) = ( J 2 )'(Coo + Coi - Cio) ■ (B.5) 
Equations (IB. 41) and (IB.5J) are identical to (IB.ll) and (1B.2I) . except with the substitution 



Cio l— ^ —Cio- In this case, we still require Coo = Coi < £2/3- This combined with ( IB. 51) 
would impose — C 10 < —2.405, where J 2 < 0.162, again guaranteeing that we will not 
surpass the optimum found in the preceding case. 



Implementing two-photon interference in the frequency domain 19 

References 

[1] Ekcrt A K 1991 Phys. Rev. Lett. 67 661-663 
[2] Waks E, Zccvi A and Yamamoto Y 2002 Phys. Rev. A 65 052310 
[3] Ma X, Fung C H F and Lo H K 2007 Phys. Rev. A 76 012307 
[4] Schcrer A, Sanders B C and Tittel W 2011 Opt. Express 19 3004-3018 

[5] Acin A, Brunncr N, Gisin N, Massar S, Pironio S and Scarani V 2007 Phys. Rev. Lett. 98 230501 
[6] Ekert A K, Rarity J G, Tapster P R and Massimo Palma G 1992 Phys. Rev. Lett. 69 1293-1295 
[7] Tittel W, Brendcl J, Zbindcn H and Gisin N 2000 Phys. Rev. Lett. 84 4737-4740 
[8] Jennewein T, Simon C, Weihs G, Weinfurter H and Zeilinger A 2000 Phys. Rev. Lett. 84 4729-4732 
[9] Naik D S, Peterson C G, White A G, Bcrglund A J and Kwiat P G 2000 Phys. Rev. Lett. 84 
4733-4736 
[10] Aspect A, Grangier P and Roger G 1982 Phys. Rev. Lett. 49 91-94 
[11] Franson J D 1989 Phys. Rev. Lett. 62 2205-2208 

[12] Tittel W, Brendcl J, Gisin N and Zbindcn H 1999 Phys. Rev. A 59 4150-4163 
[13] Ou Z Y and Mandcl L 1988 Phys. Rev. Lett. 61 54-57 
[14] Kim H, Ko J and Kim T 2003 J. Opt. Soc. Am. B 20 760-763 
[15] Fedrizzi A, Hcrbst T, Aspelmeyer M, Barbieri M, Jennewein T and Zeilinger A 2009 New Journal 

of Physics 11 103052 
[16] Shih Y H and Scrgienko A V 1994 Phys. Rev. A 50 2564-2568 

[17] Barreiro J T, Langford N K, Peters N A and Kwiat P G 2005 Phys. Rev. Lett. 95 260501 
[18] Ramelow S, Ratschbacher L, Fedrizzi A, Langford N K and Zeilinger A 2009 Phys. Rev. Lett. 103 

253601 
[19] Olislager L, Cussey J, Nguyen A T, Emplit P, Massar S, Merolla J M and Phan Huy K 2010 Phys. 

Rev. A 82 013804 
[20] Merolla J M, Mazurenko Y, Goedgebuer J P and Rhodes W T 1999 Phys. Rev. Lett. 82 1656-1659 
[21] Bloch M, McLaughlin S W, Merolla J M and Patois F 2007 Opt. Lett. 32 301-303 
[22] Capmany J and Fernandez-Pousa C R 2010 JOSA B 27 A119-A129 
[23] Capmany J and Fernandez-Pousa C R 2011 Laser and Photonics Reviews 5 750-772 
[24] Martin A, Issautier A, Herrmann H, Sohler W, Ostrowsky D, Alibart O and Tanzilli S 2010 New 

Journal of Physics 12 103005 
[25] Abramowitz M and Stegun I A 1964 Handbook of Mathematical Functions with Formulas, Graphs, 
and Mathematical Tables (New York: Dover) chap 9, p 363 ninth Dover printing, tenth GPO 
printing ed 
[26] Clauser J F and Home M A 1974 Phys. Rev. D 10 526-535 

[27] Merolla J M, Duraffourg L, Goedgebuer J P, Soujacff A, Patois F and Rhodes W T 2002 Eur. 
Phys. J. D 18 141-146 



